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RAYLEIGH  WAVES  AT  THE  CONTINENTAL  MARGIN 


by 


Julius  Kane 


abstract 


In  this  report,  the  first  of  a  series  on  the  analytic  continuation 
of  wave  functions  past  discontinuities,  we  introduce  an  elementary 
procedure  for  the  solution  of  problems  involving  the  diffraction  cf 
vector  fields.  In  particular,  the  report  discusses  the  propagation 
of  Rayleigh  waves  incident  obliquely  upon  the  continental  ma.gin.  The 
crustal  layering  on  either  side  of  the  coastline  is  modeled  mathemati¬ 
cally  as  a  two-part  boundary  layer  in  such  a  fashion  that  the 
relevant  reflection  and  transmission  coefficients  emerge  as 
elementary  algebraic  expressions.  The  procedure  that  permits  such 
a  solution  to  be  found  is  to  introduce  a  diffraction  analogue  of 
the  well-known  procedure  in  electrical  engineering:  replace  a 
transmission  line  network  by  a  lumped  parameter  equivalent  in  a 
specified  frequency  or  wave  number  domain. 
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I .  INTRODUCTION 


From  the  point  of  view  of  a  geophysicist,  mathematical 
analysis  cannot  begin  to  acquire  practical  importance  unless 
the  solution  is  in  such  a  form  that  the  answers  to  his  questions 
can  be  promptly  given.  Unfortunately,  the  class  of  elastic 
wave  diffraction  problems  which  have  tractable  solutions  usually 
describe  oversimplified  geophysical  structures.  On  the  other 
hand,  when  the  theoretician  studies  complex  problems,  too  often 
he  concentrates  on  the  mathematical  details  of  the  solution  so 
that  his  calculations  are  often  more  impressive  than  useful. 

One  powerful  mathematical  technique,  the  Wiener  Hopf  method, 
has  not  had  the  application  to  geophysics  that  it  should  have. 

The  difficulty  in  the  application  of  this  method  is  that  it 
usually  requires  an  involved  transcendental  kernel  to  be  "split" 
into  suitable  analytic  factors.  In  my  experience  with  problems 
cf  this  type,  it  has  been  my  observation  that  very  often  this 
prolix  calculation  is  of  more  mathematical  than  physical  interest^’ 
and  that  usually  the  interesting  parts  of  the  solution  depend  in  a 
very  minor  fashion  with  respect  to  this  "split"  function.  Indeed, 
the  necessity  for  this  function  -  theoretic  decomposition  can 
often  be  dropped  as  a  requirement  from  the  analysis,  and  used  as 
the  basis  of  an  approximation  procedure  to  solve  problems  beyond  the 
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scope  of  the  Wiener-Hopf  theory  .  Other  investigators  such 
as  Koiter^^  and  Carrier^  have  been  aware  of  the  relative 
unimportance  of  this  complex  factorization  and  have  suggested 
the  use  of  simpler  substitute  kernels  which  can  be  factored  by 
inspection.  While  there  is  much  merit  in  this  approach,  I  feel 
that  it  is  a  poor  philosophy  to  take  the  solution  to  an  idealized 
problem  and  then  approximate  the  answer.  A  preferable  procedure 
is  to  formulate  the  problem  in  such  a  way  that  the  mathematical 
difficulties  are  anticipated  and  avoided  from  the  start.  That 
is,  a  certain  amount  of  intuition  and  insight  into  the  nature  of 
the  problem  can  be  used  to  guarantee  that  the  subsequent  analysis 
will  lead  to  simple  representations. 


In  the  sequel,  I  should  like  to  give  an  example  of  this  philosophy 
by  considering  the  propagation  of  oblique  Rayleigh  waves  past  a 
crustal  discontinuity  such  as  the  land- sea  boundary ^^Sure  l) .  As  a  side 
condition,  an  elementary  character  for  the  solution  will  be  required. 

Of  course,  there  will  be  a  need  to  make  certain  assumptions  con¬ 
cerning  the  fields  to  be  calculated.  However,  we  shall  be  explicit 
about  these  modifications  and  indicate  how  arbitrary  improvements 
can  be  made  in  the  analysis  if  necessary. 

A  crucial  phase  in  this  study  will  be  the  need  to  formulate 
appropriate  boundary  conditions  to  characterize  the  crustal 
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layering  of  the  mantle.  In  the  simpler  problem  of  electro¬ 
magnetic  wave  propagation,  Leontovich"  ^  has  introduced  the 
notion  of  an  "impedance  boundary  condition"  to  model  the  physics 
at  an  interface.  However,  this  idea  while  very  useful  does  have 
its  limitations.  For  one  the  solution  to  the  associated  boundary 
value  problem  is  still  rather  complex,  second  the  idea  is  valid 
only  if  the  surface  impedance  is  reasonably  independent  of  the 
nature  of  the  excitation.  Although  the  impedance  concept  can 
be  generalized  to  describe  interfaces  with  more  general  properties^ 
the  intractable  nature  of  the  solution  remains.  A  feature  of  the 
subsequent  analysis  will  be  a  derivation  of  an  approximate  boundary 
condition  for  elastic  wave  diffraction  whose  merit  is  that  its 
accuracy  can  be  arbitrarily  great  without  sacrificing  simplicity 


in  the  form  of  the  solution. 
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II.  BASIC  EQUATIONS 

In  the  absence  of  body  forces,  small  displacements  S(u,v,w) 
of  an  elastic  solid  characterized  by  the  Lame  pareineters  X  ,ll  , 
and  density  p  can  be  derived  from  a  scalar  potential  $  and  a 
vector  potential 

(1)  S  =  grad  $  +  curl  ¥  (Y^O,^) 

(Three  components  in  the  vector  potential  are  redundant  for  we 
can  express  say  Y^  *  fCY^.Y^)  so  that  without  any  loss  of 
generality  we  can  set  Y^  -  0.)  For  monochromatic  vibrations, 
we  can  suppress  a  time  factor  e  ZWt  and  it  can  be  shown  that  $ 
and  the  components  Y  .  satisfy  the  reduced  wave  equations 


(V2  +  k2)  $  =  0 
c 


(V2  +  k2)  Y  -  0 

S  l 


k2  =  u)2p/(\  +  2p) 


k2  =  w2p/u 
s 


(2) 
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Simultaneous ly  we  should  like  a  solution  for  the  Rayleigh 

wave  incident  from  either  side  of  the  crustal  discontinuity. 

For  this  reason  let  v  and  v  be  the  phase  velocity  of  the 

Rayleigh  wave  on  the  left  and  right  sides  of  the  boundary,  and 

let  =  w/v^.  ar‘^  =  to/v^  designate  the  corresponding  wave 

numbers.  Also,  F  and  TD  reprerent  the  required  shear /compressional 

L  K 

ratios  for  either  side.  In  the  notation  of  figure  one,  the  variation  of 

f  101 

the  incident  Rayleigh  wave  in  the  mantle  will  have  the  fornr  J 


(3) 


r  i. 

i  me 


-  +  ay 

e 


vu 

v  inc  =  (sin  0,  0,  cos  0  ) 


rRe‘  ,llR5  +  by 


where 

(4)  ^  *  xcos  0+  z  sin  Q, 

(5)  aKkJ-k, c2)f  b.(^-kj) 


and  for  the  special  case  of  a  free  elastic  half  space  the 

coefficient  F  would  be  given  by 
K 


.  r‘  -  wH 

R  i  •  i(Vvs)2 


(6) 
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+  ^ Rz  sin  0 

The  z-variation  of  all  fields  will  be  of  the  form  c  *  ,  i . e 


(7) 


+  *  z  sin  0 

$  (x,  y .  z)  -  cp(x,  y)  e 


+  ik  z  sin  0 

Y  .(x,v,z)  *  Y  .(x,y)  e  R  ,  i  -  1,3 

z  z 


so  thatcp(x,y)  and  the  satisfy  the  reduced  two-dimensional 

wave  equations 


(8) 


djg 

d*2 

+  ±2. 

+  2 
oy 

+  p2cp  =  0 

9 

2  i2 

p  =  k 

K  c 

2  .  2  o 
kR  p \n  0 

d2Y  . 

d2'4f  . 

2 

2  ,  2 

,2  •  2  , 

1 

a,2 

V 

■r  q  >jf  .  =  o 

1  9 

q  =  k  - 
s 

kR  s  m  0 

Note  that  for  angles  of  incidence  such  that  sin  0  >  v  /v 

R  s 

2  2 

both  p  and  q  will  be  negative. 


*- plane 


fi8ure  2.  The  complex  ct  -  plane ,  The  contour  follows  the  real  axis  except  for  the  indicated 
pole  and  branch  point  deformations.  The  pole  at  l  does  not  play  any  role  in 
the  analysis  and  is  not  shown. 
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A  completely  general  solution  of  equations  (8)  can  be 
written  down  at  once  by  separation  of  variables  as 

i  r  2  2  * 

cp(x,y)  ■  2^  Jc  exP[  **x  "  *(P  -  or  )  y]  dof 

(9) 

i 

t  .(x,y)  ■  L  B(a)  exp[  *yx  -  f(q2  -  a?2)  y]  do?,  i  =  1,3 

l  4TTJ  ,Jl  t 

where  k(a)  and  B(a)  are  kernels  to  be  found,  and  the  contour  C 

i 

follows  the  real  axis  except  for  the  standard  pole  and  branch  point 
deformations  (figure  2). 

III.  THE  BOUNDARY  CONDITIONS 

Either  the  oceanic  or  the  continental  crustal  structure  can 
be  considered  to  be  a  shallow  transition  layer  that  continues 
the  field  within  the  mantle  to  the  free  boundary  at  the  earthfs 
surface.  The  details  of  the  transition  become  more  important 
for  the  shorter  period  Rayleigh  waves.  However,  the  group 
velocity  dispersion  curves  for  Rayleigh  waves  traveling  oceanic 
or  continental  paths  intersect  at  about  seventeen  seconds.  This 
implies  that  for  waves  of  period  greater  than  ten  seconds  we 
can  treat  such  a  crustal  structure  as  a  two-part  boundary  layer 
whose  coefficients  have  different  values  on  the  land  and  sea 
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sides  of  the  coastline.  For  this  purpose,  we  need  some  way  of 
relating  the  free  surface  conditions  to  the  crust-mantle  interface. 

A.  Exact  Formulation 


At  a  free  surface,  the  normal  and  tangential  stresses  must 
vanish  or 


d  u  .  d  v  .  d‘ 


yy 


dx*  oy^  oz *■  dy 


dv 


0, 


(10)  a 


xy 


du  ,  dv 

r~  +  — —  =  0, 

ay  dx 


yz 


dw  .  dv 

~  +  "T-  =  0, 

oy  dz 


which  in  terms  of  cp  and  become 

z 
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yy 


?  r^2  ^1  d 

->.kc  cp  +  2u  [-f  +  ikR  sin  e^-  -  — 5^J  -  0 

oy 


(11) 


xy 


2  2 

.  d2<p  .  .  d*l  ,  d  *3  d  *3 

■ 2  dS;  *  i\  sin  esr  +^r  •  ■° 


by  di 

a  =  zk  sin  9  [  2  +  zk  sin  0  vjr  -  — — J  -  — ~  =  0 

ayz  R  L  oy  R  1  ox  J  ^  2 


By  means  of  the  Kaskell^^  -  Thomson^^  matrix  method,  these 
equations  at  the  surface  car  be  transformed  to  an  equivalent  set  which 
refer  the  boundary  conditions  at  the  free  surface  to  the  crust- 
mantle  interface.  In  this  fashion  the  details  of  the  crustal 
layering  on  either  side  of  the  boundary  will  be  introduced  into 
the  analysis.  However,  the  boundary  conditions  so  obtained  will 
be  no  simpler  than  the  set  (11) ,  and  in  fact  a  good  deal  more 
complex.  Although  the  program  we  have  just  described  would  lead 
to  a  soluble  Wiener-Hopf  problem,  the  complexities  of  the 
analysis  would  render  this  procedure  rather  prolix  for  the 


information  desired. 
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B  The  Simplified  Formulation 

The  important  features  of  the  exact  formulation  described 
in  part  A  can  be  incorporated  in  a  much  simpler  approach.  The 
conditions  we  need  are  those  that  characterize  the  boundary 
layer  in  such  a  fashion  that  the  geometric  acoustics  poles 
mirror  the  crustal  layering  correctly.  We  c  n  assume  that  we 
..now  the  dispersion  relations  for  either  side  of  the  boundary, 
this  is  equivalent  to  making  use  of  the  Haskell  -  Thomson  matrix 
analysis  in  an  indirect  fashion. 


The  dispersion  relations  define  the  v  and  T  of  a 

K,  L  K , L 

Rayleigh  wave  uniquely.  If  we  refer  to  the  explicit  representation 
(3)  we  see  that  at  any  characteristic  depth  ^ay  y  =  0  the  Rayleigh 
wave  will  satisfy  the  set  of  equations 

dy  dy 

kR,LrR,L$  dT  '  dT  =  0 


(12) 


df 

dz  r 


R,L 
R,  L 


Y, 


=  0 


t. 

PR,L  ; 


0 
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These  equations  are  analogues  of  the  telegrapher's  equations  in 

electromagnetic  theory,  and  are  used  to  describe  the  propagation 

characteristics  of  Rayleigh  waves  in  a  layered  half  space.  They 

describe  the  surface  loading  of  the  mantle  in  an  exact  fashion 

for  Rayleigh  wave  propagation.  That  is,  they  define  the  correct 

velocity  v  ,  and  the  shear /compressional  ratio  ±  T  together 
ft >  R>  L 

with  its  proper  change  of  sign  when  the  direction  of  propagation  is 

reversed.  The  surface  parameter  f  is  depth  dependent  and 

R,  L 

transforms  according  to  the  law 


(13) 


i 

expfoc^L  -  kc2)  yDJ 
expf  (k^L  -  k2)\o] 


With  this  observation,  the  propagation  of  a  Rayleigh  wave  past  a 
step  discontinuity  can  be  described  as  a  special  case  of  the 
subsequent  analysis  if  the  step  discontinuity  is  interpreted  as 

a  jump  in  T  according  to  the  law  (13),  without  any  discontinuity 

(Figure  3) 
m  v  ° 


-15- 


r\  v 


Figure  3:  An  elastic  solid  with  a  small  step  discontinuity  can  be 
thought  of  as  a  half  space  with  different  values  of  F 
on  either  side  of  the  transition  if  the  analysis  is 
referred  to  the  reference  plane  indicated. 
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Of  course  there  are  any  number  of  other  simplified  boundary 
conditions  equations  that  the  Rayleigh  wave  will  satisfy,  but 
the  equations  in  the  set  (12)  have  several  important  advantages. 

(a)  They  are  simple  and  uncomplicated. 

(b)  They  are  isotropic  in  that  they  do  not  depend  upon 
the  Rayleigh  wave's  direction  of  propagation  in  the 
xz- plane. 

(c)  They  avoid  the  use  of  the  transverse  operator  d/dy 
which  would  have  the  effect  of  introducing  radicals 
in  the  transform  plane.  This  deliberate  avoidance  of 
the  irrational  factors  will  enormously  simplify  the 
subsequent  function- theoretic  analysis. 

(d)  They  characterize  the  crustal  layering  in  a  direct 

fashion  for  they  introduce  the  parameters  v  and 

K ,  L 

r  defined  by  the  relevant  period  equations. 

R>  L 

(e)  It  can  be  shown  that  they  imply  conservation  of  energy. 
That  is,  the  rays  which  are  reflected  from  a  boundary 
at  which  these  synthetic  conditions  are  imposed  will 
have  the  same  energy  as  the  incident  plane  wave. 
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C •  Discussion 


For  those  familiar  with  electrical  engineering  practice,  an 
analogy  between  formulations  A  and  B  may  prove  to  be  useful.  In 
simplified  terms,  the  Haskell  matrix  method  is  a  vector  analogue 
of  transforming  an  impedance  through  a  transmission  line  network 
where  the  details  of  the  crustal  structure  characterize  the 
network  parameters.  In  symbolic  terms,  the  exact  formulati.on 
can  be  considered  to  be  the  shunt  transmission  line  networks  of 
figure  3A  which  loads  a  transmission  line  representation  of  the 
mantle.  As  is  well  known,  the  impedance  characteristics  of  such 
networks  are  involved  transcendental  functions  of  frequency. 

Within  any  significant  frequency  interval  however,  a  transmission 
line  network  can  be  replaced  by  an  equivalent  network  consisting 
of  lumped  parameter  elements.  Such  networks  have  much  simpler 
impedance  functions  which  are  rational  functions  of  the  frequency. 
In  diffraction  problems,  wave  number  plays  a  role  analagous  to 
frequency,  and  the  substitute  set  of  equations  (12)  can  be 
considered  to  be  the  mathematical  equivalent  of  a  lumped  circuit 
representations  of  the  shunt  transmission  line  networks  (11)  in  a 
region  of  the  wave  number  plane  centered  about  the  Rayleigh  pole. 

From  a  mathematical  point  of  view,  the  nature  of  the  solution 
to  a  Wiener-Hopf  problem  is  that  it  analytically  continues  a  wave 
function  past  a  discontinuity.  If  two  distinct  but  reasonably 
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equivalent  sets  of  bounda  y  roi.ditions  are  used  to  characterize 
each  interface  then  the  corresponding  fields  will  also  agree 
rather  veil.  Koiter's  and  Carrier's  justification  of  the  use  of 
substitute  kernels  in  the  transform  plane  has  the  immediate 
corollary  of  validating  the  use  of  simplified  boundary  conditions 
to  describe  the  physics,  since  the  kernel  is  in  effect  the 
transform  realization  of  the  boundary  condition  operators.  Our 
use  of  the  set  (12)  as  boundary  conditions  is  however  a  preferable 
procedure  for  it  is  more  direct  and  we  can  justify  it  on  physical 
grounds  by  using  the  Haskell  method  backwards.  That  is,  we  can 
transform  the  conditions  in  (12)  from  the  crust-mantle  interface 
to  the  free  surface.  Naturally,  they  will  not  imply  that  the 
surface  stresses  vanish  except  for  Rayleigh  waves.  However,  any 
incident  ray  whose  wave  number  is  close  to  the  Rayleigh  wave 
number  will  not  induce  appreciable  stresses  at  the  surface. 

Indeed,  for  P  or  SV  waves  striking  the  free  surface  at  grazing 
angles  of  incidence,  the  equations  in  (12)  describe  the  physics 
exactly.  Essentially  all  the  diffracted  rays  are  tangential  to 
the  surface  discontinuity  and  for  waves  of  this  type  the  stress 
mismatch  at  the  surface  will  be  small.  In  any  event,  more 
terms  can  be  added  to  the  simplified  boundary  conditions  after 
the  fashion  of  the  analysis  in  [7,8,9]  to  get  any  desired  degree 
of  surface  matching. 
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IV.  SOLUTION 


In  terms  of  cp(x,y)  and  .(x,y)  the  synthetic  boundary 
conditions  (12)  at  the  crust-mantle  interface  can  be  written 


(14) 


'K?..LrR,i5’  ■  -  dT  ■  °-  y  =  0 


(15) 


iScp  +  4s 

R,  L  * 


=  0, 


V  =  0 


(16) 


x  F  3 

R,L 


=  0, 


y  «  0 


where 

(17) 


S  «  +kp  sin  a 


is  a  coefficient  that  always  involves  k  owing  to  the  suppressed 

R 

z-variation. 


Let  R(a)  denote  the  hermitian  matrix  (r^)  where  the 
coefficients  r  are  identified  as  polynomial  components  of 
the  realization  of  the  differential  operations  (14)  -  (17)  in 
the  transform  plane  when  x  is  positive,  viz 


Likewise,  let  L  (cn)  be  the  matrix  counterpart  of  (19)  but  with 
parameters  k  and  I"  appropriate  tor  the  left  boundary.  The 

Li  Li 

right  matrix  R  is  singular  when 

(19)  a  -  ±  r  =  ±  kRcos  0  , 
and  the  left  matrix  L  is  singular  when 

(20)  a  =  ±  A  =  ±  (k^  -  k“  sin^  0 ) 

These  points  correspond  to  the  geometric  acoustics  poles  for  the 
right  and  left  boundaries  respectively.  The  relationships  (19) 
and  (20)  are  seen  to  deiine  the  correct  angles  of  refraction  at 
the  interface. 


We  introduce  the  unknown  vector 
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and  a  constant  vector  m  describing  the  mismatched  components  of 
the  incident  wave  upon  the  left  boundary 


(22) 


m 


where 


(23) 


m2  “  ‘CVr/Tl  '  V  8in  0 

”3  =  !<kLr L  ’  kR)  C°8  6 
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A.  Function  -  Theoretic  Method 


Vajnsteirr  Karp1  ,  and  Clemmow1  J  have  introduced  a 

rather  general  method  for  solving  dual  integral  equations.  The 
VKC  procedure  can  be  generalized  to  solve  vector  problems  with 
either  the  exact  or  simplified  boundary  conditions. 


A  plus  (or  minus)  si perscript  on  a  matrix  or  vector  quantity 
will  indicate  that  all  of  its  components  are  analytic  in  the  upper 
(or  lower)  a- plane  above  (or  below)  the  contour  of  integration  C. 
Unless  otherwise  specified,  these  quantities  vanish  uniformly  at 
infinity  in  their  relevant  half-plane  of  regularity.  By  Jordan's 
lemma,  the  boundar>  condition  on  the  right  will  be  satisfied 


provided  that 
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(24) 


R(a)  X  =  f +(a) 


where  t  (a)  is  unknown  as  yet.  On  the  left  boundary  we  need 
satisfy 


(25) 


L(a)  X  + 


m 

a  +  r 


g'(«) 


The  term  IT)/  (a  +  k  )  represents  a  collection  of  poles  whose 
residue  contributions  represent  the  driving  effect  of  the 
incident  Rayleigh  wave  upon  the  left  boundary.  In  order  to 
solve  these  equations,  we  substitute  (24)  into  (25)  and  obtain 

(2  6)  K(cO  f +(a)  -i-  ~~ -  -  g~  (a) 

Ci  +  r  ^ 

where  the  matrix  to  be  "split"  is 

-  1 

K(a)  =  L(a)  [  R  (cr)  J 


(27) 
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For  the  simplified  boundary  conditions  the  entries  in  the  matrix 
kernel  KO2)  are  to  rational  functions  of  the  transform  variable 
Ct  so  that  by  elementary  moans  we  can  split  K(a)  as 

_1  + 

(28)  K (0!)  =  [N  (a)]  p  (a) 

For  the  exact  formulation,  however,  contour  integration  would  bo 
required  to  perform  the  factorization.  In  any  event,  the  matrix  factors 
P+(0!)  and  N  (d)  are  regular  in  their  respective  half  planes,  however  un 
lixe  the  vectorsf +(q;)  or  9  ( PC, )  they  do  not  vanish  at  infinity 
but  approach  a  constant  value. 

If  we  add  and  subtract  the  same  term,  we  can  rearrange 
equation  (26)  into  the  form 


(29)  P+(a)  f +(a)  + 


N'j-r) 


a  4-  r 


=  N"  (pc)  g"  (a)  + 


a  +  r 


[  N  (-  - )  -  N  (a) 1  ~n 


Each  side  cf  this  equation  is  regular  in  a  half  plane  and  defines 
the  analytic  continuation  of  an  entire  function.  Owing  to  the 
assumed  growths  of  f+(3)  and  9  (cO  this  entire  function  vanishes 
uniformly  at  infinity  and  must  thus  be  the  null  constant.  Hence 
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we  can  set  each  side  of  (29)  to  zero  and  obtain 

,  -1  - 

[P  fa)]  N  (-r>m 

f  (a)  - - - 


(30) 

-1 

IN  fa)] 

9  fa)  -  -  ----- —  [N  <-r)  -N  fa)]  rn 


With  the  aid  of  the  first  of  these  expressions,  the  desired 
solution  X  is  seen  to  be 


[Rfa)]  [P+fa>]  N’(-r) 


(31) 


X 


a  +  r 
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B.  Alf abraic  Method 

For  rational  kernels  the  general  procedure  described  in 
Section  A  is  rather  abstruse  and  a  more  direct  approach  can 
yield  the  solution  in  a  simple  algebraic  fashion;  at  the  same 
time  the  transparency  of  the  analysis  will  offer  an  insight 
into  the  mathematics  that  will  suggest  a  generalization  that 
can  be  used  to  solve  tnnre  complex  problems. 

The  matrices  L,  and  R  are  such  that  the  only  singularities 
we  can  have  in  the  solution  X  are  poles  in  the  complex  plane  at 
a  ~  ±  r ,  ±  JL.  We  know  that  the  vector  X  must  have  entries  which 
are  rational  functions  of  a.  The  poles  of  X  can  only  arise  at 
a  *=  ±  r  (incident  and  reflected  waves),  and  at  a  **  -i  (transmitted 
wave)  for  the  residue  wave  at  d  =  4-  Si  is  a  non-physical  contri¬ 
bution.  Hence,  X  must  be  inversely  proportional  to  (a  -  r)  (a  +  r) 
(pt  +  H)  ,  and  we  can  immediately  concentrate  attention  on  the 
numerators  of  its  components.  That  is^  X  must  be  of  the  form 
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(32) 


X 


_ l _ 

(a  -  r)  (a  +  r)  (a  +  i) 


I  »!  +  >2(a  '  r)  \ 

*TR  sin  0  [a1  +  a3(a  -  r)  (a  -  a^/a^)  ] 

l  , 

\  ffR  cos  e  [  a1  +  a5  (OC  -  r)  ]  f 


where  the  are  five  constants  to  be  determined.  The  order  of 
the  polynomial  entries  in  X  is  easily  established  if  we  anticipate 
that  we  need  use  Jordan's  lemma  after  multiplication  by  either 
the  operator  R  or  1.  .  The  five  unknown  constants  a  .  can  and  must  be  so 
chosen  that  all  boundary  conditions  are  obeyed.  For  x  positive, 
we  need  satisfy 


(33) 


1 

2rr  i 


KXe  ***  d a 

Je 


o, 


x  > 


0 


and  the  form  of  (32)  is  such  that  it  automatically  meets  this 
constraint,  since  the  only  residue  contribution  in  the  upper 
half  plane  is  a  reflected  i  eigh  wave. 
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For  the  boundary  condition  on  the  left  side  we  need  satisfy 


(34) 


■At  J<i 


m  .  fox  . 

)  e  do;  =  0, 


a  4-  r 


x  <  0 


and  by  Jordan's  lemma  this  requirement  will  be  et  if  the 
residue  coefficients  at  CL  =  -r,  -l  are  null.  That  is  we  need 
satisfy 

(35) .  (ct  +  r)  L(-r)  X(-r)  +  m  *  0, 

and 

(36)  (a  +  jO  LC-i)  X(-£)  =  0, 

Equation  (36)  is  a  homogeneous  relation  and  can  be  satisfied  by 
specifying  two  constants  to  make  the  residue  contribution 
proportional  to  a  transmitted  Rayleigh  wave  that  matches  the 
left-hand  conditions,  viz. 


sin  H  [  a^  +  (l  +  r)  (j l  +  a^/a^)]  =-fL  sin  [a^  -  a^i-t  +  r)  ] 


(37) 


Tr  cos  e  [ax  -  a_  (A  +  r)  ] 


=+fL  cos  0t  [at  -  (l  +  r)] 
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The  three  equations  in  (35)  represent  the  left  boundary's  reaction 
that  cancels  the  incident  fields  and  together  with  the  two  in  (37) 

we  determine  the  five  constants  a  .  uniquely. 

i 
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V.  DISCUSSION  OF  THE  SOLUTION 


The  elementary  character  of  the  solution  (31)  or  (36) 
permits  a  ready  evaluation  of  the  scattered  fields.  For 
y  =  0  and  x  positive  the  entire  field  is  a  collection  of 
residue  contributions  which  can  be  grouped  into  an  incident 
Rayleigh  wave,  and  a  reflected  Rayleigh  wave  with  a  relative 
amplitude 


(38) 


ft  = 


2r  (r  +  £) 


Our  convention  will  be  to  describe  the  coefficient  of  the 
compressional  potential  of  the  Rayleigh  wave  as  its  amplitude. 
Although  our  solution  includes  scattered  body  waves,  their 
contribution  is  absent  for  y  =  0  owing  to  our  formulation. 


On  the  left  sid*3  of  the  boundary,  the  incident  Rayleigh  wave 
is  cancelled  by  pole  contributions  as  it  should.  The  remaining 
residue  terms  are  seen  to  represent  a  transmitted  Rayleigh  wave 
of  amplitude 


a1  -  a2(r  +  &) 
2r (£  -  r) 


(39) 
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This  wave  travels  in  a  direction  Q  which  is  given  by 

VL 

»in  Q  = - sin  0  , 

VR 

a  familiar  expression  for  the  refraction  of  a  ray  at  an  inter¬ 
face.  In  general  both  phase  velocities  v  and  v  depend  upon 
frequency,  and  as  a  consequence,  the  angle  0^  is  frequency 
dependent.  Thus  the  diffraction  of  a  Rayleigh  pulse  incident 
obliquely  upon  a  crustal  interface  will  sort  the  transmitted 
spectral  components  into  different  directions  after  the  fashion 
of  an  optical  prism. 
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